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ABSTRACT 

2002cx-like supernovae are a sub-class of sub-luminous Type la supernovae. Their light 
curves and spectra are characterized by distinct features that indicate strong mixing 
of the explosion ejecta. Pure turbulent deflagrations have been shown to produce such 
mixed ejecta. Here, we present hydrodynamics, nucleosynthesis and radiative transfer 
calculations for a 3D full-star deflagration of a Chandrasekhar-mass white dwarf. Our 
model is able to reproduce the characteristic observational features of SN 2005hk (a 
proto- typical 2002cx-like supernova), not only in the optical, but also in the near- 
infrared. For that purpose we present, for the first time, five near-infrared spectra of 
SN 2005hk from —0.2 to 26.6 days with respect to i?-band maximum. Since our model 
burns only small parts of the initial white dwarf, it fails to completely unbind the 
white dwarf and leaves behind a bound remnant of ^ 1.03 Mq - consisting mainly of 
unburned carbon and oxygen, but also enriched by some amount of intermediate-mass 
and iron-group elements from the explosion products that fall back on the remnant. 
We discuss possibilities for detecting this bound remnant and how it might influence 
the late-time observables of 2002cx-like SNe. 

Key words: supernovae: individual: SN 2005hk - techniques: spectroscopic ~ meth- 
ods: numerical - hydrodynamics - radiative transfer 



1 INTRODUCTION 

Due to an empirical relation between their peak bright- 
ness and their light curve evolution ( |Phillips|1993j Pskovskii| 
19771, Type la supernovae (SNe la) can be used as stan- 
dardiseable candles to measure the expansion history of the 
Universe (|Riess et al.|p98l [Schmidt et "al]p98} [Perlmut- 



ter et al. 19991. However, we still have no detailed picture 



about the progenitor systems of these luminous explosions 
beyond the general idea that they originate from the ther- 
monuclear disruption of carbon-oxygen (CO) white dwarfs 
(WDs) (for a review see e.g. [Hillebrandt fc Niemeyer|2000[ ). 



* Based on observations collected at ESO, Parana! Program ID 
076.D-0183 

f E-mail; mkromer@mpa-garching.mpg.de 



Radioactive isotopes, in particular ^®Ni and its daughter nu- 
cleus ^^Co, which are synthesised within the thermonuclear 
burning, power the observable display of SNe la. 

As of today several sub-classes of SNe la have been iden- 
tified (Li et al. 20111. Among the most peculiar of these 



events are the SN 2002cx-like objects (Li et al. 20031. For 



their light curve decline rate, these objects are sub-luminous 
with respect to the Phillips relation and do not show sec- 
ondary maxima in the near-infrared (NIR) bands. Their 
spectra are characterised by very low expansion velocities 
compared to other SNe la and show signs of strong mixing 
in the ejecta ( |Jha et al.|[2006[ [Phillips et aL]|2007[ ). This is 
in clear contrast to normal SNe la, which are characterised 
by strongly layered ejecta ( Stehle et al.|20(j5 Mazzali et al. 
|2007) . 

While explosion models involving a detonation are not 
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able to explain such an ejecta structure (e.g. Sim et al 
2010 Seitenzahl et al.|2012 1, such a strong mixing can easily 



be obtained from turbulent deflagrations in WDs ( Gamezo, 



Khokhlov, & Oran 2004 Ropke & Hillebrandt 20051. The 



low expansion velocities are also in good agreement with 
the small amount of kinetic energy releeised in deflagration 



models. Such considerations led Branch et al. (20041 and 



Jha et al. (20061 to conclude that SN 2002cx-like objects 



might be related to pure deflagrations of Chandrasekhar- 
mass WDs. This interpretation was supported by |Phillips| 



et al. (20071, who compared the broad-band light curves of 



SN 2005hk (a well-sampled proto-typical 02cx-like event) to 
synthetic light curves of a 3D deflagration model of |Blin-| 



nikov et al. ( 2006 1 and found good agreement between the 



data and the model. 

The models presented by [Blinnikov et al.| ( |2006| l , how- 
ever, do not allow for a detailed comparison to the observed 
spectral time sequence since the multi-group approximation 
employed in their radiative transfer simulations is too coarse. 
Moreover, their underlying hydrodynamic explosion models 
are restricted to one spatial octant of the progenitor WD, in- 
troducing artiflcial symmetries to the flame evolution. Here, 
we report on a 3D full-star deflagration simulation through 
to the homologous expansion phase. We perform detailed ra- 
diative transfer calculations with the time-dependent 3D ra- 
diative transfer code ARTIS ( Kromer fc Sim|200"9 Sim|2007 1 
from which we obtain a time series of synthetic spectra that 
we compare to the observed spectra of SN 2005hk. 

The paper is organised as follows. In Section|2]we give a 
brief description of our explosion simulation and present the 
resulting ejecta structure. In Section[3]we present synthetic 
observables for this explosion model and compare them to 
the observed light curves and spectra of SN 2005hk. Finally, 
we discuss our results and give conclusions in Sections [4] and 
[5] respectively. 



2 EXPLOSION SIMULATION 

A Chandrasekhar-mass WD is believed to undergo about 
a century of convective carbon burning in the centre be- 
fore a thermonuclear runaway finally sets in. Since this so- 
called simmering phase is characterised by highly turbu- 
lent flows, it cannot be fully accounted for by present-day 
numerical simulations (but see e.g. 'Hoflich & Stein 



Kuhlen, Woosley, & Glatzmaier 2006; Zingale et al. 



2002 



2009 



Nonaka et al. 12012 1. Thus, the actual ignition configuration 



of Chandrasekhar-mass WDs is not well constrained. Given 
this ignorance, one may use the ignition geometry as a free 
parameter and explore a larger set of explosion simulations 
with various ignition setups. Thereby one has to account for 
both different ignition strengths and ignition positions. A 
good way to achieve the former is to use a multi-spot ig- 
nition scheme which seeds unstable burning modes in a ro- 
bust and numerically well-controlled way. Recently, we have 
performed such a systematic study for different ignition se- 
tups of 3D full-star pure deflagration simulations (Fink et 
al., in preparation) yielding ^^Ni masses between 0.035 and 
0.38 M©. 

Here, we focus on a detailed comparison to SN 2005hk. 
For that purpose we select one of the models of the series 
by Fink et al. that produces a ^''Ni mass on the order of 




300 km 



Figure 1. Ignition setup of the N5dcf model. Shown is a vol- 
ume rendering of the WD in mint green colour. As discussed in 
Section [2] the ignition kernels were randomly placed within a ra- 
dius of 150 km around the centre of the WD, as illustrated by the 
bluish sphere. The exact configuration of the ignition kernels is 
shown in the enlarged inlay (see also Table [TJ. 



the observationally derived value of SN 2005hk. Applying 
Arnett 's law (|Arnett|198"2 1 to the observed bolometric light 
curve, iPhillips et al.l (120071) report a '^'^Ni mass of ~O.2M0 



for this SN. With a '^^i mass of 0.18 Mq, model N5def of 
the Fink et al. series comes close to that value. 

In the N5def simulation, an isothermal (T = 5 x 10^ K) 
Chandrasekhar-mass WD was set up in hydrostatic equi- 
librium with a central density of 2.9 x 10^ gem"'' and a 
homogeneous composition of carbon and oxygen in equal 
parts by mas^ To account for an assumed solar metallic- 
ity of the zero-age main-sequence progenitor, we start with 
a of 0.49886, corresponding to 2.5 per cent of ^^Ne in 
the initial composition. The WD was then discretized on a 
three-dimensional Cartesian moving grid (R6pke'2005) with 
512^ cells consisting of two nested parts (central resolution 
1.9 km) and ignited in flve spherical ignition kernels that 
were placed randomly in a Gaussian distribution within a ra- 
dius of 150 km from the WD's centre. By chance, for model 
N5def this algorithm produced a fairly one-sided ignition 
configuration with all kernels lying in a relatively small solid 
angle. Thus, this configuration is representative for a slightly 
off-centred single-spot ignition with a larger number of ini- 
tially excited burning modes. The actual configuration is 
shown in Fig. [l] All kernels have a radius of rk — 10 km and 
are at distances d between 6.3 and 80.9 km from the origin 
(see Table [1]). 

Neglecting any possible defiagration-to-detonation 
transition (see [Seitenzahl et al.|2012 for an alternative evo- 
lution of this model in a delayed-detonation scenario) we fol- 



^ Although the detailed flame evolution depends on the exact 
value of the central density (which is not well constrained for 
Chandrasekhar-mass WDs, e.g. Seitenzahl, Ciaraldi-Schoolmann,| 
|fc Ropke„20TTt , it is not expected to have a qualitative impact 
on the outcome of the explosion. 



© 2012 RAS, MNRAS 000,[T|fT2l 




2002cx-like SNe from deflagrations with bound remnants 3 




Figure 2. Snapshots of the hydrodynamic evolution of our model N5def. Shown are volume renderings of the mean atomic number 
calculated from the reduced set of species in the hydrodynamic simulation (see colour bar). To allow a view to the central part of the 
ejecta, a wedge was carved out from the front of the ejecta. (i) At 0.75 s a one-sided deflagration plume rises towards the WD surface 
and fragments due to Rayleigh- Taylor and Kelvin-Helmholtz instabilities, (ii) At 1.5 s the expansion of the WD quenches the burning 
and the explosion ashes wrap around the unburned core, (iii) Finally, at 100 s the unburned core is completely engulfed by the explosion 
ashes which are ejected into space. The small triads at the bottom right corner of each panel indicate the scaling at the origin of each 
plot; from left to right the legs of the triads represent 1000 km, 4000 km and 500 000 km, respectively. 



Table 1. Position of the ignition kernels of model N5def. Given 
are the x, y and z coordinates of the centre of the individual 
ignition kernels and their distance d to the centre of the WD. 



# 


X 


y 

(in 


z 

km) 


d 


1 


65.5 


-15.5 


24.0 


71.5 


2 


38.6 


-22.7 


67.3 


80.9 


3 


13.0 


8.2 


15.1 


21.6 


4 


-5.0 


-51.3 


-2.6 


51.7 


5 


5.6 


2.9 


0.6 


6.3 



lowed the flame evolution up to 100 s after ignition when the 
ejecta reach homologous expansion. To this end we used the 
LEAFS code, a three-dimensional finite-volume discretiza- 
tion of the reactive Euler equations which is based on the 



PRO METHEUS implementation ( [Fryxell, Miiller, fc Arnett 
1989[) of the 'p iecewise parabolic method' (PPM) by lColella 



& Woodward ( 1984 1. Deflagration fronts are modelled as dis- 
continuities between carbon-oxygen fuel and nuclear ash, 
and their propagation is tracked with a level-set scheme 



( Reinecke et al.||1999 Osher fc Sethian||1988 Smiljanovski, 



Moser, fc Klein|1997 \ . All material crossed by these fronts is 
converted to nuclear ash with a composition and energy re- 
lease depending on fuel density. Composition and energy re- 
lease are interpolated from tables which have been calibrated 
using our full 384-isotope network. After a very short phase 
of laminar burning following ignition, the propagation of 
deflagrations is dominated by buoyancy- and shear-induced 
instabilities and interactions with a complex turbulent flow 
field. The unresolved acceleration of the flame due to turbu- 
lence is accounted for by a sub-grid scale model (Schmidt,' 
Niemeyer, fc Hillebrandt|2006a[ [Schmidt et al\2006hy Self- 



The hydrodynamic evolution of our model is shown 
in Fig. [2] Since a deflagration flame cannot burn against 
the density gradient, our asymmetric ignition conflguration 
leads to the formation of a one-sided deflagration plume that 
fragments due to Rayleigh- Taylor and Kelvin-Helmholtz in- 
stabilities. Once the deflagration front reaches the outer lay- 
ers of the WD, the burning quenches due to the expansion of 
the WD and the ashes wrap around the still unburned core 
until they finally engulf it completely. A similar evolution of 
the deflagration flame was already described for single-spot 



off-centre ignitions by e.g. Plewa, Calder, fc Lamb (20041 
and [Ropke, Woosley, fc Hillebrandt| ( |2007| ). While |Plewa| 
et al. ( 2004 1 found an ensuing detonation to be triggered 



when the ashes collide on the far side of the star (see also 
|Seitenzahl et al.|2009al), w e - similarly to |Ropke et al.| ( |2007] ) 
and Jordan et al. (2012b I - do not find high enough densi- 



ties and temperatures for such a detonation to occur due to 
a significant expansion of the WD during the deflagration 
phase. 

Since only a moderate fraction of the core of the WD 
is burned, the nuclear energy release in our simulation 
{Enuc ~ 4.9 X 10^" erg) is less than the binding energy of 
the initial WD (Sbind = 5.2 x 10^" erg). Nevertheless, about 
0.37 M© of the WD are accelerated to escape velocity and 
ejected into the ambient medium with a kinetic energy of 
1.34 X 10^'^ erg. The remainder of the mass of the initial WD 
is left behind and forms a bound remnant. Similar findings 



gravity is dealt with by a monopole gravity solver. 



were already reported in 2D (e.g. Livne, Asida, fc Hofiich 
2005 1 and recently in the context of a failed gravitationally- 
confined detonation ( Jordan et al.|2012b K 

To obtain detailed nucleosynthesis yields of the explo- 
sion, we performed a post-processing calculation with a 384 
isotope network for 10® Lagrangian tracer particles which 
were passively advected during the hydrodynamic simula- 
tion to record thermodynamic trajectories of mass elements 
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Table 2. Yields of select species for model N5def. 





Bound remiiciiit 


Ejecta 


Total 


1.028 


0.372 


C 


0.422 


0.043 


o 


0.484 


0.060 


Ne 


0.054 


0.005 


Mg 


0.004 


0.013 


Si 


0.015 


0.025 


S 


0.004 


0.009 


Ca 


0.0003 


0.001 


Fe 


0.004 


0.031 


Ni 


0.025 


0.187 


56 Ni 


0.022 


0.158 



hydro dynamic simulations (by which point homologous ex- 
pansion is a good approximation) to a 50^ Cartesian grid 
and followed the propagation of 10* photon packets for 111 
logarithmically-spaced time steps between 2 and 120 d after 
explosion. To speed up the initial phase, a grey approxima- 
tion, as discussed by Kromer & Sim (20091, was used in 



optically thick cells, and the initial 10 time steps (i.e. the 
initial three days after the explosion) were treated in local 
thermodynamic equilibrium (LTE). For our simulation we 
used the 'big_,5f-4' atomic data set of |Kromer fc Sirii (20091 
with a total of ~ 8.2 x 10® bound-bound transitions. The 
resulting broad-band light curves and a spectral time series 
are shown in Figs. [4] and [5] and compared to SN 2005hk as 
a proxy for SN 2002cx-like objects. 



( Travaglio et al.|2004 Seitenzahl et al.|2010 1. A compilation 
of the masses of important species is given in Table|2] To de- 
termine the mass that stays bound in the remnant and that 
of the unbound ejecta, we calculated the asymptotic specific 
kinetic energy tkin.a = ekin,f + tgrav.f for all tracer particles. 
Here, ekin.f = Uf/2 and tgrav.f are the specific kinetic and 
gravitational binding energies sXt = 100 s, i.e. at the end of 
our simulation, respectively. Only for positive values of ekin.a 
a particle will be able to escape the gravitational potential. 
Otherwise, it will stay bound. 

Finally, we mapped the unbound tracer particles to 
a Cartesian grid to reconstruct the chemical composition 
of the explosion ejecta in the asymptotic velocity space[^ 
To this end, we used an SPH like algorithm as described 
in [Kromer et alT] ( |2010[ ). The resulting ejecta structure is 
shown in Fig. |3] In contrast to the simulation by |Jordan| 
al. (2012b I, our model shows no pronounced large-scale 
composition asymmetry. There are, however, small scale 
anisotropies due to the turbulent evolution of the defiagra- 
tion fiame. The density jump at about 2000 km s~^ is pro- 
duced during the strongest pulsation of the bound core (at 
t ~ 12 s): After significant initial expansion due to burn- 
ing, the mostly unburned inner parts of the ejecta start 
falling back inwards. This leads to the formation of a central 
compact core, which, after maximum compression, starts to 
expand again. At the same time, outer layers are still in- 
falling and an accretion shock forms at the edge of the dense 



3.1 Broad- band light curves 



core (cf. Bravo & Garcfa-Senz 20091. In our model, maxi- 
mum temperatures in the accretion shock just reach 10^ K 
at densities of 5 x 10® gem"'', therefore no detonation will 



be triggered (for critical conditions see Ropke et al. 2007 
and [Seitenzahl et al.||2009b[ ) . 



3 SYNTHETIC OBSERVABLES 

To obtain synthetic spectra and light curves for our model, 
we used the time-dependent 3D Monte Carlo radiative trans- 
fer code ARTIS ( [Kromer fc Sim|2009| [Sim|2007) . For the ra- 
diative transfer simulation we mapped the abundance and 
density structure of the unbound ejecta at the end of the 



For unbound tracer particles, the asymptotic velocity is deter- 
mined as = ,y2€^~. 



As reported by Jha et al. (20061 and Phillips et al. (20071, 



SN 2002cx-like objects are characterised by distinct spectral 
features and Kght curves. In the following, we will investi- 
gate to which extent our model is capable to reproduce these 
features. From the broad-band fight curves in Fig. |4| it is 
immediately obvious that our model matches the low peak 
luminosities of 2002cx-like SNe quite well. This, of course, is 
not so surprising, since we picked a model with a ^''Ni mass 
close to that observationally derived for SN 2005hk. How- 
ever, we do not only get the correct peak luminosities but 
also approximately the right colours at maximum. Moreover, 
our model naturally explains the absence of secondary max- 
ima tn the NIR bands. This is a consequence of the turbulent 
burning, which leads to an almost completely mixed ejecta 
structure with a roughly constant iron-group element (IGE) 
mass fraction at all velocities (see Figs. |2] and |3|. [Kasen 



(20061 has shown that SN ejecta with such a homogenized 



composition are not expected to show secondary maxima in 
the NIR. 

Other peculiarities of 2002cx-like SNe are their slowly 
declining light curves in R and redder bands and conse- 
quently also in UVOIR bolometric. In this respect our model 
has some shortcomings. While it nicely reproduces the post- 
maximum decline in the U, B and V bands, our post- 
maximum light curve evolution is too fast particularly in R 
and redder bands, but also in UVOIR bolometric. This could 
be a consequence of our low ejecta mass of only 0.37 Mq, 
which is not able to trap enough 7-radiation after maximum. 
We note, however, that the synthetic fight curves in Fig. [4| 
are only powered by radionuclides in the ejecta. Thereby 
we neglect any possible influence from the bound remnant, 
a puffed-up stellar object heated by the explosion, and also 
from the ®®Ni that falls back onto the remnant (see Table[2|. 
After maximum fight, when the ejecta start to become opti- 
cally thin, the bound remnant will be exposed and it could 
kick in as an additional luminosity source. We will discuss 
this intriguing possibility in Section [4. 3| 

However, looking at the overall shape of the early fight 
curves, an insufficient trapping of radiation in the ejecta 
and thus too low ejecta mass seems to be more likely. With 
a B-band rise time of 11.2 d, our model evolves somewhat 
too fast compared to the rise time of SN 2005hk, for which 
( PhUlips et aL][2007t find a value of 15 ± 1 d. 
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Figure 3. Composition of tlie unbound ejecta of model N5def. Sfiown are slices through the x-z-plane in the asymptotic velocity space 
at 100 s after the explosion. The top left panel displays the density, the following panels show the mass fraction of selected species (C, 
O, Si, stable iron-group elements and ^^Ni from top left to bottom right). 




-10 10 20 30 40 -10 1 20 30 40 -10 1 20 30 40 -10 1 20 30 40 

Days relative to B 



Figure 4. Synthetic light curves of our model. The panels from top left to bottom right contain UBVRIJHK bolometric and broad-band 
U ,B ,V ,R,I ,J,H light curves ( [Bessell fc Brett|1988[ |Bessell|1990^ . The dark red line corresponds to the angle-average of the model. The 
light red region indicates the scatter caused by 100 different viewing angles. Time is given relative to B-band maximum. The black 
upside-down triangles show the observed light curves of SN 2005hk ( [Phillips et "al]|2007| . Due to missing /f-band data for SN 2005hk 
the bolometric light curve shown only accounts for UBVRIJH fluxes. For comparison, the grey symbols show observational data of three 
well-observ ed normal SNe la: SN 2 003du (circles, [Stanishev et al.|2007| , SN 2004eo (squares, [Pastorello et al.|2007a[ , and SN 2005cf 
(triangles, |Pastorello et al.|2007bl l. 
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4000 5000 6000 7000 8000 9000 



Rest wavelength in A 

Figure 5. Spectral evolution of model N5def from -6.3 to 14.6 
days (epochs are given with respect to B-band maximum). For 
comparison the observed spectra of the 2002cx-like SN 2005hk 
( [Phillips et al.|2007| l are over-plotted in black. The observations 
were de-reddened and de-redshifted according to the values of 
|Phillips et"al] | |2007| l and |Sahu et al.| | |2008| l, respectively. 



3.2 Optical spectra 

|Jha et al.| ( |2006[ ) and [Phillips et"aL] ( |2007[ ) also report pe- 
culiar features for the spectra of 2002cx-like SNe. At pre- 
maximum epochs they find an SN 1991T-like spectrum with 
a blue continuum, weak absorption features at ~4200 and 
5000 A and no strong Call absorption features. As can be 
seen from the top panels of Fig. [5j our model naturally 
predicts this behaviour, characteristic of an ionisation state 
higher than in normal SNe la. 

Moreover, 2002cx-like SNe show a significant contribu- 
tion from IGEs to the spectra at all epochs, while normal 
SNe la do not show strong IGE features until well-after 
maximum light. Due to the strong mixing of the ejecta, 
our model is almost homogenized and thus predicts strong 
IGE contributions at all epochs. This is illustrated in Fig.|6] 
There, the colour coding below the synthetic spectrum in- 
dicates the fraction of escaping packets in each wavelength 
bin that were last emitted by bound-bound transitions of 
a particular element (the associated atomic numbers are il- 
lustrated in the colour bar). From this it is clearly visible 
that emission by IGEs already dominates the spectrum be- 




2000 5000 10000 20000 

Rest Wavelength (A) 



Atomic number 



2 6 7 8 9 10 1 1 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 

Figure 6. Formation of the UV to NIR SED for our N5def model 
for three different snapshots. Epochs are given relative to B-band 
maximum. The colour coding indicates the elements responsible 
for both bound-bound emission and absorption of quanta in the 
Monte Carlo simulation. The region below the synthetic spec- 
trum is colour-coded to indicate the fraction of escaping quanta 
in each wavelength bin that were last emitted by a particular ele- 
ment (the associated atomic numbers are illustrated in the colour 
bar). Similarly, the coloured regions along the top of the plots in- 
dicate which elements were last responsible for removing quanta 
from a particular wavelength bin (either by absorption or scatter- 
ing / fluorescence) . White regions between the emerging spectrum 
and the colour-coded bound-bound emission indicate the contri- 
bution of continuum processes (bound-free, free-free) to the last 
emission. 



fore maximum light, while this is not the case for stratified 
models, such as e.g. W7 (see figure 2 of [Kromer fc Sim|2009[ ). 

Another characteristic feature of the spectra of 2002cx- 
like SNe are extremely low line velocities at all epochs, indi- 
cating significantly lower expansion velocities and thus less 
kinetic energy than in normal SNe la. Since our model burns 
only parts of the WD, the kinetic energy is rather low com- 
pared to other explosion models that fully burn and unbind 
the progenitor WD. As a consequence, the line velocities of 
our model are in good agreement with the observed ones 
(compare Fig. [5|. 
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1.0x10' 1.4x10' 1.8x10' 2.2x10' 

Rest wavelength in A 

Figure 7. NIR spectra of the 2002cx-like SN 2005hk (black) and 
synthetic spectra of model N5def (red) from —1.3 to 26.6 days 
(epochs are given with respect to S-band maximum). The ob- 
servations were de-reddened and de-redshifted according to the 
values of [Phillips et al.| l |2007| l and |Sahu et alT] | |2008| l, respec- 
tively. For comparison, NIR spectra of the normal SN la 2005cf 
| |Gall et al.|2012] l are shown for similar epochs (blue). 



3.3 Near-infrared spectra 

In Fig.[7]we show, for the first time, NIR spectra of a 2002cx- 
like SN. These spectra of SN 2005hk were collected by the 
European Research Training Network (RTN) 'The Physics 
of Type la Supernova Explosions]^ (details on the observa- 
tions and data reduction are given in Appendix [A|). Com- 
pared to normal SNe la there are a few prominent differ- 
ences. While in normal SNe la the characteristic Co ll emis- 
sion humps in the K band (e.g. [Gall et al.||2012"| do not 
form before 20 d after maximum light, they emerge at 6.6 d 
after maximum in SN 2005hk (see grey highlighted region in 



http://www.mpa-garching.mpg.de/ ~ rtn/ 



Fig.[7|. At later epochs (~ 30 d) the Coll features appear sig- 
nificantly more narrow in SN 2005hk than in normal SNe la. 
A similar but even stronger effect can be observed in the J 
and H bands, where very narrow spectral features develop 
at these epochs. 

Our model naturally reproduces these characteristic 
features of the NIR spectra of SN 2005hk. The narrow fea- 
tures at later epochs are a direct consequence of the low 
kinetic energy released during the explosion and reflect the 
low expansion velocities of the eject a (Fig. [3]) compared to 



that of models of normal SNe la (e.g. Nomoto, Thielemann, 
fc Yokoi|[l984l [Ropke et ar]|2012| ). The early occurrence of 



the Con emission humps can be understood as a result of 
the almost completely mixed ejecta. While IGEs are con- 
centrated in a central core in models of normal SNe la, the 
turbulent burning of our deflagration model leads to a very 
homogeneous ICE distribution in the ejecta, which extends 



to the highest velocities. As shown by Kasen (20061, such a 
mixing does not only lead to single peaked NIR light curves 
but also to an increased NIR flux around B-band maximum. 

Although our synthetic NIR spectra reproduce the char- 
acteristic features of SN 2005hk, a detailed comparison 
shows some differences. While the absolute flux level and the 
overall SED agree quite well with the observations around 
maximum light, differences start to emerge about a week 
after maximum light. At these epochs the flux level of our 
model is too low compared to the observations as already 
discussed for the broad-band light curves (Section |3.1[ ). 



4 DISCUSSION 

4.1 Spectropolarimetry 

There are a few observational constraints which we cannot 
address directly with our synthetic observables. |Chornockl 
et al. (20061 and Maund et al. (20101 presented spectropo- 



larimetric observations of SN 2005hk and reported a line and 
continuum polarization on a level of a few tenths of a per- 
cent, which is typical of SNe la ( Wang fc Wheeler|2008 1 and 
requires approximately spherical explosion ejecta. Since our 
radiative transfer simulations currently do not include the 
polarization state of the radiation fleld, we cannot directly 
compare the observed polarization spectra to our model. We 
note, however, that our ejecta structure does not show very 
pronounced large-scale asymmetries (see Fig. [3| and that 
the synthetic observables are not very sensitive to different 
lines-of-sight as visible from the light curves in Fig. |4] Thus, 
a weak polarization signal seems plausible. 



4.2 Late-time spectra 

Further observational constraints come from late-time (~ 9 
months after explosion) spectra, where 2002cx-like SNe show 
extremely narrow permitted Fell lines and a continuum or 



continuum like SED (Jha et al. 2006 Phillips et al 



20071. In contrast, normal SNe la are in the nebular phase 



at these epochs. Jha et al. ( 2006 1 also report a tentative 



detection of weak Ol features in the late-time spectra of 
SN 2002cx. They do, however, not see strong features of [O l] 
AA6300,6364 which have been predicted for 3D deflagration 



models (Kozma et al. 20051, where unburned O is present 



© 2012 RAS, MNRAS 000,(1^121 



8 M. Kromer et al. 



down to the lowest ejecta velocities due to turbulent burning 
as in our model. 

Since ARTIS does not yet account for non-thermal ion- 
ization and excitation, which are the dominant processes 
at late epochs (e.g. Kozma & Fransson 1998a|b I, we can- 
not directly test the late-time spectra of our models against 
the observations. Whether [O l] features will arise, however, 
depends strongly on a possible microscopic mixing of dif- 
ferent species. This cannot be addressed with present-day 
numerical models, since it takes place on scales which are 
not resolved. However, if such a mixing is present, transi- 
tions of other species stronger than [O l] AA6300,6364 will 
dominate the cooling. Also the central density of our explo- 
sion ejecta (~120gcm~^ at 100 s after explosion) is about 
30 times larger than in the W7 model 



' 4 g cm 



Nomoto 



et al. 19841 or typical delayed detonation models that are 



in good agreement with normal SNe la. This might be im- 
portant in explaining the peculiar features of the late-time 
spectra of 2002cx-Uke SNe. 

Beyond that, we can only speculate that the bound rem- 
nant in our model might be able to explain the peculiar non- 
nebular features of the late-time spectra of 2002cx-like SNe 
if the observed spectrum was a superposition of emission 
from the bound remnant and the SN ejecta. We will discuss 
the implications of the bound remnant in more detail in the 
next section. Alternatively, Sahu et al. ( 2008 1 tried to ex- 



plain the observed late-time spectra of SN 2005hk by com- 
bining synthetic photospheric (e.g. Mazzali & Lucy 1993 1 
and nebular (e.g. Mazzali et al.|p001 1 spectra of the ejecta 



of a parametrized ID explosion model. 



4.3 Implications of the bound remnant 

As discussed in Section[2] our model does not burn the com- 
plete WD but leaves behind a bound remnant. Due to the 
strong expansion of the ejecta, at late times this remnant 
cannot be resolved from our simulations, thus preventing 
any detailed analysis. However, it is clear that the remnant 
must be a puffed- up stellar object heated during the explo- 
sion. From the bound tracer particles we can derive the over- 
all composition of the remnant (see Table |2| , which shows 
that it is enriched with IGE material and some amount of 
radioactive ^^Ni. This enrichment is due to fall-back of ex- 
plosion ashes onto the bound remnant, so that the enriched 
material should be located in the outer layers of the remnant 
(see also |Jordan et al.|2012b|| . 

We also note that the ""Ni content of the bound rem- 
nant is large enough to contribute to the optical display of 
SN 2005hk at late epochs (Fig. [8| if 7-rays are completely 
trapped in the remnant. Since the remnant does not expand, 
its densities are significantly larger than those of the ejecta, 
making this assumption plausible. Having some contribution 
from the bound remnant at late epochs, could also explain 
the highly peculiar spectra of 2002cx-like SNe, which never 
become fully nebular (see discussion in previous section). In 
particular the extremely narrow IGE features might be ex- 
plained if the emission originated from an IGE-rich crust of 
the bound remnant. 

Being heated by the actual explosion and long-lived ra- 
dioactive isotopes, the bound remnant might even become 
directly detectable at much later epochs when the explosion 
ejecta go into the supernova remnant phase. For SN 2005hk, 




100 200 300 

Time since explosion in days 



400 



Figure 8. Bolometric light curve of SN 2005hk (black diamonds; 
following j Phillips et al.] l |2007^ , we adopted a rise time of 15 d). 
For comparison, the red line shows the synthetic bolometric light 
curve from our radiative transfer simulation in the ejecta. At 
t = 120 d our detailed non-grey radiative transfer simulation ends. 
Thereafter, we extended the synthetic light curve with a calcu- 
lation using a grey UVOIR opacity. The blue curve shows the 
instantaneous energy deposition due to the ^^Ni decay sequence 
in the bound remnant, assuming full 7-ray trapping. The magenta 
curve shows the sum of both these contributions. 



observations were performed with the Gran Telescopio Ca- 
narias on July 20, 2012, i.e. 2443.5 d after _B-band maximum. 
These, however, have not shown any point source down to a 
limiting magnitude of Mr — —8.08 ± 0.49 (for details of the 
observation see Appendix [B| . 

After the bound remnant has radiated away all the heat 
from the explosion and possible reheating from radioactive 
isotopes, it will again become a WD, however with a peculiar 
composition enriched in IGEs (and also intermediate-mass 
elements). Moreover, [Jordan et al. ( 2012b| l report kick ve- 
locities of up to 520kms~ for the remnant WDs in their 
failed gravitationally-confined detonation simulations, thus 
claiming the existence of hyper-velocity enriched WDs as 
a 'smoking gun' for this scenario. From our simulations we 
do not see such kicks. This difference may originate from 
the different gravity solvers used. In general, no perfect mo- 
mentum conservation is to be expected from approximate 
solutions of the Poisson equation. 



4.4 Implications on the binary companion star 

One problem of single-degenerate progenitor models dis- 
cussed in the literature is that of the non-detection of hydro- 
gen lines in the observed (nebular) spectra of SNe la (e.g. 
Leonard|2007 1. Recent hydrodynamical simulations investi- 
gating the impact of SN la ejecta on a non-degenerate com- 
panion star (e.g. Pan, Richer, Taam|2012 Liu et al.|2012 \ 
have shown that the amount of hydrogen stripped from the 
companion star seems to be significantly larger than the ob- 
servationally derived upper limits (Leonard 20071. Due to 



the low kinetic energy of our explosion model, compared 
to that of normal SNe la, only a small amount of material 
should be stripped from the companion in the progenitor 
systems of SN 2002cx-like SNe, thus probably avoiding a 
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signature of hydrogen lines in late-time spectra. This should 
be addressed by detailed impact studies in the future. 

4.5 The class of 2002cx-like SNe 



After the discovery of SN 2002cx by |Li et al.| (|2003|) as a 
peculiar object, it soon became evident that this was not a 
unique event but rather the prototype of a new class of pecu- 



liar SNe la with a striking spectral homogeneity ( Jha et al 



|2006[ ). From a volume-limited sample of the Lick Observa- 
tory Supernova Search (LOSS), Li et al.| ( |2011 1 estimate that 
SN 2002cx-like explosions contribute at about 5 per cent to 
the total SN la rate. 

Recently, SN 2008h a ( [Foley et al.|2009| ) and SN 2007qd 
I McClelland et al.|2010 | were proposed as additional mem- 
bers of the class. Though spectroscopically similar to 
SN 2002cx, those objects are fainter by more than 2 magni- 
tudes than a typical 2002cx-like SN and show far lower line 
velocities, on the order of 2000 km s~^ , while SNe 2002cx and 
2005hk showed line velocities on the order of 7000 km s^^. 



Nevertheless, McClelland et al. (20101 found a relationship 



between light-curve stretch, the peak brightness and the ex- 
pansion velocities among their (small) sample of objects. 
This indicates that their extended class of SN 2002x-like 
objects may originate from a single explosion mechanism. 

With an observationally derived ^^Ni mass of 0.003 Mq 
for SN 2008ha ( [Foley et al.|2009[ ), however, this object can- 
not easily be explained in our model of a deflagration of 
a Chandrasekhar-mass WD. From a systematic study of 
3D full-star pure deflagration simulations for different ig- 
nition setups (Fink et al., in preparation) we obtain a mini- 
mum ^®Ni mass of 0.035 M©, which is more than a factor 10 
larger than the observationally derived value for SN 2008ha. 
Though we did not sample all possible ignition setups in the 
Fink et al. study nor take into account different progenitor 
compositions, this difference may be difficult to reconcile 
with a deflagration of a Chandrasekhar-mass WD. Thus an- 
other explosion mechanism might be at work for these faint 



objects. Foley et al. (2009 20101 suggested deflagrations of 



sub-Chandrasekhar-mass WDs, while Valenti et al. (20091 



and Moriya et al. (20101 favoured a core-collapse origin. 



4.6 Rates 



As mentioned in the previous section, Li et al. (20111 es- 
timate from the LOSS that 2002cx-like SNe contribute at 
about 5 per cent to the total SN la rate. This is remarkably 
close to some theoretical predictions for the contribution 
of single-degenerate hydrogen-accreting systems to the total 
SN la rate. From binary population synthesis models, e.g. 
Ruiter, Belczynski, & Fryer (2009) find a value of ~3 per 
cent (their model 1 with acE x A = 1). At the same time 
their predicted total SN la rate for all progenitor channels 
considered (single-degenerate Chandrasekhar-mass scenario 
and double-degenerate mergers with a total mass above the 
Chandrasekhar-limit) is comparable to the Galactic SN la 
rate (Ruiter et al.|2009, figure 2). Taking this at face value, 
it could be that all single-degenerate systems will end up as 
2002cx-like SNe. Normal SNe la would then originate from 



other progenitors like doublc-dcgoncrate mergers (e.g. Pak 



mor et al.[[2012 Ruiter et al. 


j2012) or sub-Chandrasekhar- 


mass double detonations (e.g. 


Fink et al.|2010 Kromer et al. 



2010 Woosley fc Kasen|2011 1. However, the observationally 
derived SN rates as well as the binary population synthesis 
rate predictions for the single-degenerate scenario are some- 



what uncertain. For example, Han & Podsiadlowski ( 2004 1 



can reconcile the observed Galactic SN la rate with the birth 
rate of SNe la in the single-degenerate Chandrasekhar-mass 
scenario from calculations with their binary population syn- 
thesis code which assumes a different prescription for hydro- 



gen accretion than Ruiter et al. ( 2009 1 



4.7 The fate of Chandrasekhar-mass explosions 

In a companion study (Fink et al., in prep.) we are currently 
investigating the outcome of a larger set of 3D deflagration 
models. Within this sample we find a wide range of explo- 
sion strengths depending on the actual ignition conditions. 
However, to obtain ^^Ni masses on the order derived for 
SN 2002cx or SN 2005hk, only a very narrow range of igni- 
tion configurations works. 

We find that appropriate deflagrations with moderate 
^^Ni production result from ignition setups that can be in- 
terpreted in the scheme of a slightly off-centre single-spot 
ignition. Thus, if all Chandrasekhar-mass explosions led to 
2002cx-like SNe, as hypothesized in the previous section, 
only such ignition configurations should be realized in Na- 
ture. Note that this is in agreement with state-of-the-art pre- 
ignition simulations (e.g. Kuhlen et al.[[2006 Zingale et al 
2009[[Nonaka et al.|20"l2| ). 

However, it could also be that for other configura- 
tions an ensuing detonation after the initial defiagration 
phase occurs, thereby producing a layered ejecta structure 
as needed for normal SNe la (e.g. [Mazzali et al.[[2007[ ). 
Several mechanisms to trigger such a secondary detona- 
tion have been proposed: a spontaneous deflagration-to- 
detonation transitio n ([Khokhlov[ [1991[ ), a gravitationally 
confined detonation ( Plewa et al.|2004[), a pulsating reverse 
detonation (Bravo & Garcfa-Senz''2006) or a pulsationally- 
assisted gravitationally-confined detonation ( [Jordan et al.[ 
2012a I. However, for none of these models the formation of 



the secondary detonation can be resolved in detail. Instead, 
detonations are triggered when certain conditions are met 
which are thought to be sufficient to ignite a detonation 
from simplified ID simulations. 



5 CONCLUSION 

We have presented hydrodynamic simulations and radiative- 
transfer calculations for a 3D pure-defiagration of a 
Chandrasekhar-mass WD. Our particular model (N5def), 
chosen from a larger sample of simulations (Fink et al., in 
preparation) , has an asymmetric ignition configuration lead- 
ing to a relatively weak deflagration, which fails to com- 
pletely unbind the WD. Only ~O.37M0 of material are 
ejected, whereas the remainder stays gravitationally bound, 
leaving behind a stellar remnant of ~ 1 Mq . 

The ejecta contain a sig nificant '^^Ni mass (~0.16Mq), 
have low kinetic energy, are well-mixed and relatively sym- 
metric on large scales. All of these properties make N5def a 
promising model for 2002cx-like SNe la. Moreover, a com- 
parison of our synthetic observables with SN 2005hk shows 
good agreement. Our synthetic light curves reproduce the 
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peak luminosity and colour of SN 2005hk and naturally pre- 
dict the lack of a secondary maximum in the NIR bands. 
Our synthetic spectra are also in good agreement with ob- 
servations. At optical wavelengths, N5def shows the same 
high ionisation, low expansion velocities and early domi- 
nance of IGE lines observed in SN 2005hk. In addition, we 
have presented and compared with previously unpublished 
NIR spectra of SN 2005hk, allowing us to verify that the 
good match between model and observations extends into 
the NIR regime. The main shortcoming of the model is an 
overly rapid post-maximum fading at red and NIR wave- 
lengths that is also apparent in the bolometric light curve, 
which is too narrow. This points at too low opacities, which 
could be the result of insufficient ejected mass. 

The most intriguing feature of the model presented here 
is the bound remnant. Unfortunately, at late times the rem- 
nant is not resolved in the simulation, so that its exact state 
after the explosion remains largely unknown. We do find, 
however, that some fallback of nucleosynthetically processed 
material (including ~ 0.02 M© of ^^Ni) onto the remnant 
will occur, resulting in a peculiar IGE-rich composition of 
the remnant object. Heated by the explosion, this will be 
a puffed-up stellar object. We can also speculate that the 
bound remnant may be sufficiently luminous to contribute 
to the observed display of the SN, explaining in part the 
highly peculiar late-time spectra of 2002cx-like SNe, which 
do not become fully nebular even one year after the explo- 
sion. 

While N5def seems to be a reasonable model for 2002cx- 
hke SNe with '^^Ni masses of ~O.2M0, it appears unlikely 
to produce ^^Ni masses two orders of magnitude lower in the 
deflagration of a Chandrasekhar-mass WD. These would be 
needed to explain 2008ha-like SNe, which have been sug- 
gested to be fainter and less energetic siblings of 2002cx-like 
SNe, originating from the same explosion mechanism. 

Finally, we point out that the observed rate of 2002cx- 
like SNe (about 5 per cent of the total SN la rate) is 
quite similar to the contribution of the single-degenerate 
Chandrasekhar-mass explosion channel to the total SN la 
rate as expected from some binary population synthesis cal- 
culations. It is thus tempting to speculate that deflagration- 
to-detonation transitions might actually not happen in 
SNe la, and that single-degenerate Chandrasekhar-mass ex- 
plosions might only lead to pure deflagrations with a 2002cx- 
like optical display. 



APPENDIX A: OBSERVATION LOG AND 
DATA REDUCTION OF NEAR-INFRARED 
SPECTRA 

Five low-resolution NIR spectra of SN 2005hk were obtained 
at VLT/ISAAC (Table|AT]|. During each night when the SN 
was observed, one of the following bright stars Hip 1115, 
Hip 26093 or Hip 115347, all of spectral type B3V, were also 
observed close in time and airmass to the SN observations. 
Four different instrument settings were used to obtain wave- 
length coverage between 0.97 mm and 2.5 mm. The observa- 
tions for each setting were reduced separately. The spec- 
tra from the different instrument settings that did not over- 
lap were scaled so that the synthetic J and H photometry 



Table Al. Log of the VLT/ISAAC NIR spectroscopy 



Date 


JD 


Phase 


(UT) 




(day) 


2005/11/09 


2 453 683.58 


-1.4 


2005/11/14 


2 453 688.58 


3.6 


2005/11/17 


2 453 691.58 


6.6 


2005/11/23 


2 453 697.58 


12.6 


2005/12/07 


2 453 711.58 


26.6 



The observations were performed in ABBA sequences, 
where A and B denote two different positions along the slit. 
The images were bias and fiat field corrected, and trimmed. 
Because of the field distortion of ISAAC, spectra taken at 
two different slit positions are not parallel. Besides, the im- 
age of the slit is also curved. For the subsequent data reduc- 
tion it was necessary to correct for these two effects. This 
was achieved with the TWODSPEC.LONGSLIT package in 
IRAF. The input data for deriving this correction were: (i) 
the arc-lamp spectra obtained during the nights when the 
SN was observed and (ii) special (STARTRACE) observa- 
tions obtained by ESO twice per year and consisting of spec- 
tra of a bright star taken at 11 different positions along the 
slit. After applying the geometric distortion correction the 
night sky lines were parallel to the image columns and the 
spectra to the image rows. Next, for each pair of AB images 
the B image was subtracted from the A image. The negative 
spectrum was shifted to the position of the positive one and 
subtracted from it, resulting in an image with the sum of 
the spectra but without the sky background. All such im- 
ages were summed into a single image and the ID spectra 
were then optimally extracted. We note that the optimal ex- 
traction algorithm has to be applied on images with the pixel 
levels in the form of actual detected counts, and so it does 
not work quite correctly if applied to background-subtracted 
images. Special care was thus taken to calculate the opti- 
mal extraction weights correctly. The wavelength calibration 
was done with arc-lamp spectra. Its accuracy was checked 
against the sky lines and small corrections were applied if 
necessary. 

To remove the strong telluric absorption features, the 
SN spectra were divided by the high signal-to-noise spec- 
tra of the B3V stars. To obtain the relative flux calibra- 
tion, the result was multiplied by a model continuum spec- 
trum for a B3V star. We used a model spectrum with 
left = 18500 K and logg = 4 by R. Kurucz, available at 
|http://kurucz. harvard.edu/ The residual features due to 
the hydrogen absorption lines of the B3V stars were removed 
by dividing the SN spectra by the composite spectrum for 



the B3V spectral type from Pickles (19981. The composite 



B3V spectrum was first normalized to the continuum and 
smoothed to the instrumental resolution of each of the four 
instrumental setups used. 



APPENDIX B: GTC LATE-PHASE 
PHOTOMETRY 



matched the observed one published by Phillips et al. ( 2007 1 



Deep r-band (Fukugita e t al.|1996| ) images of the SN 2005hk 
explosion site with exposure times of 3 x 200 s were ob- 
tained at the Gran Telescopio Canarias (GTC) equipped 
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with OSIRIS on UT 2012 July 20, 2443.5 d after B-band 
maximum. The images were debiased, flat-field corrected, as- 
trometrically aligned and averaged. An image quality of 0.9 
arcsec was measured from the full width at half maximum of 
isolated stars. The photometric zero point of the image was 
determined by comparing instrumental magnitudes of stel- 
lar objects with the apparent r-band magnitudes reported 



in the SDSS data release 6 catalogue ( Adelman-McCarthy 



et al. 20081. Since no source was visible at the position of 



SN 2005hk, we estimated the limiting magnitude in our im- 
age by measuring faint stellar sources in similar locations 
as SN 2Q05hk. The apparent magnitude of the faintest of 
these sources, r = 25.68±0.41, was then adopted as Hmiting 
magnitude for the remnant of SN 2Q05hk. Note that this is 
~ 3.9 mag fainter than the last (Bessell) _R-band magnitude 



of SN 2005hk, taken 381.6 d after _B-band maximum (|Sahu 
2008[). With a colour excess of E{B 



et al. 



0.112 

mag and a distance modulus fi — 33.46 ±0.27 mag, our Hm- 
iting apparent magnitude translates into a hmiting absolute 
magnitude of = -8.08 ± 0.49. 
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